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Abstract 

We study indirect CPT violating effects in meson decays and mixing, taking 
into account the recent constraints on the CPT violating parameters from the Belle 
collaboration. The life time difference of the Bj meson mass eigenstates, expected 
to be negligible in the standard model and many of its CPT conserving extensions, 
could be sizeable (~ a few percent of the total width) due to breakdown of this fun- 
damental symmetry. The time evolution of the direct CP violating asymmetries in 
one amplitude dominated processes (inclusive semileptonic B^ decays, in particular) 
turn out to be particularly sensitive to this effect. 
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The suggestion for two distinct lifetimes for the Bd or B s meson mass eigenstates 
originated in parton model calculations [|IJ, which, at that time, were limited by numer- 
ous uncertainties of hadronic the bag parameter, top quark mass, ...) and weak 
parameters (CKM matrix elements). Many of these, however, cancel in the ratio 

where Am^Ar,;) is the mass (width) difference of the Bd meson mass eigenstates, r] t , rj 
are calculable perturbative QCD corrections, x t = and 

f f x ) = lnx-(- + - - -r-^] • (2) 

JK J 2(l-x) 3 V 4 4,(1 -x) 2(l-x) 2 J KJ 

Following the discovery of mixing in the Bd system |2|], Am^ was measured and m t 
was the only major source of uncertainty in the ratio. Using the then lower bound on m< 
it was shown |3|, |j] that AI^ is indeed very small, while Ar s , the width difference of the 
B s meson mass eigenstates could be rather large as is indicated by the scaling law j3j 
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where Vi/s are the elements of the CKM matrix and 



X Bq = (B q \[q^{l- 1 ,)bf\B q ). (4) 



In the meanwhile, many advances have taken place with the discovery of the top 
quark and the determination of its mass || and more precise values for CKM matrix 
elements. Combining the new values with the above scaling laws, the width difference 
among the Bd states is (AT/T)d « 0.0012, which is unobservable, but for B s eigenstates 
(Ar/r) s ~ 0.045. More recent calculations using heavy quark effective theory and im- 
proved QCD corrections || |7j suggest that calculations based on the absorptive parts of 
the box diagram improved by QCD corrections give reasonable estimates for both Bd and 
B s systems. 
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Nevertheless the possibility that there are loopholes in the above calculations cannot 
be totally excluded. For example, Ar g (q = d or s) is determined by only those channels 
which are accessible to both B q and B q decays. Its computation in the parton model may 
not be as reliable as the calculation of Y q , the total width which depends on fully inclusive 
decays and quark- hadron duality is valid. 

In addition to the expected phenomena, one should, therefore, be prepared for unex- 
pected effects and the final verdict on this subject should wait for experimental deter- 
mination of Ar from the B-factories, B-TeV or LHC-B. Many different suggestions for 
measuring Ar s have been put forward |3|, §. It is believed that (AY/Y)d ~ 0.1 can be 
measured at B-factories B while (AY/Y) d ~ 0.001 pi might be accessible at the LHC. 



In this article we wish to emphasize that apart from dynamical surprises in the de- 
cay mechanism, a possible breakdown of the CPT symmetry contributes to AY/Y. The 



currently available constraints on CPT violating parameters |TT] , [T2|| certainly allow this 
possibility. If this happens its effect will be more visible and detectable in the B d system 
which, in the electroweak theory, is expected to have negligible {AY /Y) d - In other words 
the scenario with {AY /Y) d large not only due to hitherto unknown dynamics but also 
due to a breakdown of CPT is quite an open possibility. In the case of (Ar/r) s CPT 
violation may act in tandem with the already known electroweak dynamics to produce an 
even larger effect. 



There are several motivations for drawing out a strategy to test CPT symmetry. From 
the experimental point of view all symmetries of nature must be scrutinized as accurately 
as possible, irrespective of the prevailing theoretical prejudices. It may be recalled that 
before the discovery of CP violation, there was very little theoretical argument in its 
favour. 

There are purely theoretical motivations as well. First of all the CPT theorem is valid 
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for local, renormalizable field theories with well defined asymptotic states. It is quite 
possible that the theory we are dealing with is an effective theory and involving small 
nonlocal/ nonrenormalizable interactions. Further the concept of asymptotic states is 
not unambiguous in the presence of confined quarks and gluons. It has been suggested 
that physics at the string scale may indeed induce nonlocal interactions in the effective 



low energy theory leading to CPT violation [TJ1. Moreover, modification of quantum 



mechanics due to gravity may also lead to a breakdown of CPT [14 



One of the major goals of the B-factories running at KEK or SLAC is to reveal CP 
violation in the B system. The discrete symmetry CPT has not yet been adequately 
tested for the B meson system, although there are many interesting suggestions to test 



it [|T|, [16]. In all such works, however, the correlation between Ar and CPT violation 
was either ignored or not adequately emphasized. It will be shown below that Ar can in 
general be numerically significant even if CPT violation is not too large. 

We consider the time development of neutral mesons M° (which can be K° or D° or B^ 
or Bg) and their antiparticles M°. The time development is determined by the effective 
Hamiltonian Hij = — with and being the dispersive and absorptive parts 
of the Hamiltonian, respectively |17| . CPT invariance relates the diagonal elements 



Mil = M 22 and r u = T 22 . (5) 

A measure of CPT violation is, therefore, given by the parameter 

^ _ H 22 - Hn 
\/H 12 H 21 

which is phase convention independent. In order to keep the discussion simple we shall 
study the consequences of indirect CPT violation only. Since indirect CPT violation 
is a cumulative effect involving summations over many amplitudes, it is likely that its 
magnitude would be much larger than that of direct violation in a single decay amplitude. 
It is further assumed that CPT violation does not affect the off-diagonal elements of H^. 



These assumptions can be justified in specific string models [I3|, where terms involving 
both flavour and CPT violations receive negligible corrections due to string scale physics. 
A further consequence of this assumption is that the usual SM inequality Myi Ti 2 
holds even in the presence of CPT violation. 



The eigenf unctions of the Hamiltonian are defined as 



\M X ) =p 1 \M°) +q 1 \M°) and \M 2 ) = p 2 \M°) - q 2 \M°) 
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with the normalization |pi| 2 + \qi\ 2 = \p2\ 2 + |<7 2 | 2 = 1- We summarize the consequences 
of the symmetries. We define 



Qi 

Vi = — 

Pi 



Q2 

m = — 

P2 



1+ 4. 



s 2 V' 2 s 



+ 2 



H- 



21 



1+ 4. 



5 2 \ 1/2 5 



H\2. 

H21 



H 



12 



1/2 



1/2 



(9) 



and note that CPT violation is contained in the first factor, while indirect CP violation 



is in the second factor with the square root. In many expressions we need the ratio 
10 — 771/772 = which is only a 
Im u = 0, Re us = 1 and 771 = rj 2 . 



10 — 7/1/772 = which is only a CPT violating quantity. CPT conservation requires 



The time development of the states is determined by the eigenvalues 
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which can be parametrized as A12 = 2 — \ Ti, 2- The quantities which occur in the 
asymmetries are \\ — A 2 = Am — |Ar and T = + r 2 ). To leading order in (T12/M12) 
they are expressed in terms of the CPT parameter 

/ an 1/2 
y= 1+4 



as follows: 

1 r 12 

Am = mi — m2 = 2|M 12 |(Re y H — Re Im y) 

2 M\2 

Ar = r 1 -r 2 = 2|M 12 |(Re -^Re y - 21m y) . (11) 

iWl2 

In the CPT conserving limit y = 1 and the contribution to Am is large, overwhelming 
CPT violating corrections. The CPT conserving contribution to Ar, on the other hand, is 
suppressed by Re j^. The purely CPT violating term dominating Ar remains, therefore, 
an open possibility. In order to get a feeling for the magnitude of AT/T, we use the small 
5 approximation and obtain | AT/r| = 0.5 x (Am/T) x (Re 5 x Im S). 

Most of the measurements of Am/T have been carried out by assuming CPT conser- 
vation. If CPT is violated its magnitude could be somewhat different (see, e.g., Kobayashi 
and Sanda in [|15j). Recently the Belle collaboration has determined Am with and without 
assuming CPT symmetry |12|| . The two results , Am = 0.463 ±0.016 and 0.461±0.008 
± 0.016 ps^ 1 respectively, do not differ appreciably from each other or from the aver- 
age Am given by the particle data group(PDG). We shall, therefore, use throughout the 
paper Am/T = 0.73, which is perfectly consistent with the PDG value. The relevant 
limits on CPT violating parameters from Belle are Im^o — m^ol/m^o <1.6 xl0~ 14 and 
|r#o — r^ |/r_B0 < 0.161, which implies |Re 5\ < 0.54 and |Im S\ < 0.23. A choice like 
Re 8 x Im 5 ~ 0.1, consistent with the above bounds, would then yield AT/T of the order 
of a few %, larger than the SM estimate by an order of magnitude. Moreover AT/T will 
be well within the measurable limits of LHC B, should 5 happen to be much smaller. 

The Belle limits are derived under the assumption that AT/T is negligible. We em- 
phasize that for a refined analysis of CPT violation Am/T and AT/T along with 5 should 
be fitted directly from the data. Such a combined fit may open up the possibility that 
5 could be somewhat larger than the above bounds. In our numerical analysis values 
consistent with the bounds as well as somewhat larger values will be considered. 
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The time development of the states involves, now, the time factors 

/-(*) 



/+(*) 
/+(*) 



-i\\t 



+ toe 



+ e 



-i\ 2 t 



and 



(12) 
(13) 
(14) 



The new feature is the presence of the factor u in the second and third of these equations. 

The decays of an original \B°) or \B°) state to a flavor eigenstate |/) vary with time 
and are given by 



P f (t) = <f\B°(t)> =\f + (t)\ 2 N <f\B°> 



P r (t) = < f\B°(t) > =/+(*) N < f\B° > 



P f (t) = < f\B°(t) > = M 2 \Mt)\ 2 N < f\B° > 



P f (t) = < f\B°(t) > =\f4t)\ 2 N\uj\ 2 <f\B°> /\ m 



(15) 
(16) 
(17) 



where iV -1 = |1 + uo\ 2 and the matrix elements on the right-hand side g = (f\B°), g = 
(f\B°), . . ., are computed at t — and have no time dependence. From these expressions 
it is evident that the five unknowns, T, Am, Ar and Re 5 and Im 5 (or equivalently 
Re uj and Im u), must be determined from the data. We emphasize that Ar must be 
treated as a free parameter, since in addition to the CPT violating contributions it may 
also receive contributions from new dynamics. In addition taking linear combinations of 
these decays, we can produce exponential decays accompanied by oscillatory terms which 
help in separating the various contributions. It may be recalled that the time dependent 
techniques for extracting these probabilities and the associated electroweak parameters 
from data are now being used extensively. 



Different schemes for testing CPT violation suggested in the literature |I5|, |TjJ often 
involve observables specifically constructed for this purpose. Here we wish to point out 
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that some of the observables involving the above probabilities, which are now being rou- 
tinely measured at BABAR and BELLE are also sufficiently sensitive to CPT violation 
and have the potential of either revealing the breakdown of this fundamental symmetry 
or improving the limit on the CPT violating parameter. One such observable is the direct 
CP violating asymmetry in Bd and Bd decays to flavor specific channels / and /, respec- 
tively fl8 |, but with / different from /. The following ratio is at the center of current 
interest 
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(19) 



In the SM or in any of its CPT conserving extensions, /+(£) = f+(t) and the asymmetry 
is time independent in general. The time independence holds even if Ar happens to be 
large due to new dynamics or direct CPT violation and/ or new physics influence the 
hadronic matrix elements. Time evolution of this asymmetry is, therefore, a sure signal 
of indirect CPT violation. Flavour specific B decays involving a single lepton or a kaon 
in the final state are possible candidates for this measurement. 



This consequence is even more dramatic for decays dominated by a single amplitude in 
the SM, in which case \g\ = \g\ and a^(t) vanishes at all times. Purely tree level decays 
arising from the subprocess b — > UiUjdk (i ^ j), penguin induced processes b — > didid k , 
dominated by a single Penguin operator or inclusive semileptonic decays b — ► Xl + u (1 = 
eor/i and X is any hadronic final state) are examples of such decays. The last process is 
particularly promising. A single amplitude strongly dominates the decay not only in the 
SM but also in many extensions of it. The large branching ratio (~ 20% for 1 = e and fi ) 
and reasonably large efficiency of detecting leptons is sufficient to ensure the measurement 
of this asymmetry at B - factories, provided it is of the order of a few percent. 
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For this class of decays the matrix elements along with their theoretical uncertainties 
cancel out in the ratio. Consequently in presence of indirect CPT violation the time 
dependent asymmetry is the same for all one-amplitude dominated processes and the 
statistics may be improved by including several channels. If a difference in the time 
dependence of various modes is observed, the assumption of one amplitude dominance 
will be questionable and new physics beyond the standard model leading to \g\ ^ \g\, in 
addition to indirect CPT violation, may be revealed. 

In Figure 1 we present the asymmetry for a one amplitude dominated process as a 
function of time for Im 5 = 0.1 and Re 5 = 0.1 (solid curve, here AT/T = 0.004) , 
0.5 (dotted curve, AT/T = 0.02) and 1.0 (dashed curve, Ar/r=0.04). Both the time 
evolution and the nonvanishing of the asymmetry are clearly demonstrated. 

The correlation between a^p and Ar calls for a more detailed analysis. As has been 
noted Ar is significantly different from the SM prediction only if Im 8 x Re 5 ^ 0. The 
numerator and the denominator of the asymmetry are determined to be 

D(t) = Pf(t)-Pf(t) 

~(M 2 - l) (e- V2t - e- Vlt ) - 4 Im uj e~ rt sin Ami] N (20) 



and 



S(t) = Pf(t) + P f (t) 



uj\ 



|2 



+ l) (e- V2t + e~ Vlt ) + 4 Re we" ri cos Amt] N (21) 



A non- vanishing asymmetry can arise in various ways 

i) Im uj t^O, which requires Im<5 7^0.0, 

ii) Rew^ 1.0 and AT/T as small as in the SM, 

or from a combination of the two possibilities. It is trivial to express uj in terms of 5 
and confirm that both D(t) and S(t) are modified from the SM prediction through 8. 
When both numerator and denominator of the asymmetry are measured accurately, one 



can determine separately real and imaginary parts of delta. This may indicate, albeit 
indirectly, that Ar is unexpectedly large. 



In order to have an idea of how large the effects can be, in Figure 2 we plot D(t) as a 
function of time for the values Re 5 = Im 5 = 0.1 (the solid curve). D(t) vanishes for 
Im 5 = and has a relatively weak dependence on Re 5, as illustrated by also plotting on 
the same figure the cases with Re 5 = 0.5 (the dotted curve) and 1.0 (the dashed curve). 
A similar study of S(t) is presented in Figure 3. This quantity is fairly insensitive to Im 
5. 

In order to estimate roughly the number of tagged B-mesons needed to establish a 
non-zero D(t) , we assume that at t=0 there is a sample of Nq tagged B® and B®. Let of 
number of semileptonic B® (B® ) decays in the time interval t=(1.0 ± 0.1) x tb be n(t) 
{n(t)) (we assume the lepton detection efficiency to be ~ 1). By requiring 

\n(t) -n(t)\ 



> 3.0, (22) 



we obtain for Re 5 = 0.1 and Im 5 = 0.1, iVo ~ 2.0 x 10 6 , a number which is realizable at 
B - factories after several years of run and certainly at the LHC Including other flavour 
specific channels like B® — > K + + X, which has a larger branching ratio (~ 70%), a 
measurable asymmetry may be obtained with a smaller N . 

In the presence of indirect CPT violation, the time integrated asymmetry is obtained 
by integrating the numerator D(t) and the denominator S(t). This leads to 

4 Re u \ 



//i ,2 ,\ Ar 4 Im uj x\ I //. l2 _ \ 4 



(23) 



with x = Am / T. In the standard model and for processes dominated by one amplitude 
the integrated asymmetry vanishes. In extensions of the SM in which the decays are no 
longer dominated by one amplitude the integrated asymmetry may be nonzero [[OJ . Thus 
a nonzero integrated asymmetry points either to new physics (coming from additional 
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amplitudes) or to indirect CPT violation. In Figure 4 we present the variation of this 
observable with Im 5 for Re 5 = 0.1 (solid line) and 0.75 (dashed line). 



Experimental studies of CPT violating phenomena can be combined with experiments 
that search for a AT/T. For example, one can consider untagged B mesons decaying to a 
specific flavour ^ |J. The observable 

S&) = P f (t) + P f (t) 

which in the absence of CPT violation have a time dependence governed by two exponen- 
tials. If now CPT violation is also included, then an oscillation is superimposed on the 
exponentials. The original articles [|3], f| considered B s decays but the same properties 
hold for Bd meson decaying semileptonically or to specific flavour final states. 

Looking at flavour non-specific channels there are results for a^ P (also denoted by 



C^) from BABAR @ and BELLE for the channel B tt+tt". In the SM using 



naive factorization this asymmetry turns out to be small f22|| . It is interesting to note 
that although the Babar result is fairly consistent with the SM prediction, the BELLE 
result indicates a much larger asymmetry. It should, however, be noted that there are 
many theoretical uncertainties. Neither the magnitude of the penguin pollution nor the 
magnitude of the strong phase difference between the interfering amplitudes can be com- 
puted in a full proof way. Direct CP violation in flavour specific, charmless decays have 
also been measured [[23|. Here the data is not yet very precise and the theoretical uncer- 
tainties are also large. In view of these uncertainties it is difficult to draw any conclusion 
regarding new physics effects. This underlines the importance of inclusive semileptonic 
decays which are theoretically clean and the branching ratios are much larger than any 
of the above exclusive modes. 

Bd decays to CP eigenstates have been observed and established CP-violation via 
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time dependent measurements p4"| , The golden example is B° —>■ ipK s where the time 
dependent asymmetry is proportional to sin 2/3, where (3 (also denoted by 0i) is an angle 
of the unitarity triangle. The current averaged value of this parameter is sin 2/? =0.78 ± 
0.08. . It is straight-forward to obtain the asymmetry in the presence of indirect CPT 
violation. An attempt to fit the data as in the SM would lead to an effective sin 2/3 which 
is time dependent. We have checked that with Re 5 = 0.1 and Im 5= 0.1, this effective 
sin 2/3 varies between 0.74 and 0.84. We therefore conclude that if sin 2/3 is determined 
with an accuracy of 5 % or better, some hint of indirect CPT violation may be obtained. 
However, this observation cannot establish CPT violation unambiguously. Since CPT 
conserving new physics may change the phase of the — Bd mixing amplitude and/or 
the decay amplitudes and lead to similar effects. 

Many other observables specifically constructed for the measurement of CPT violation 
|T5| , |I6| have been suggested in the literature. It will be interesting to compare the 
sensitivities of these observables to the CPT parameter S with that of the observables 
considered in this paper, which are already being measured in the context of CP violation. 

In summary, we wish to emphasize again that an unexpectedly large life time differ- 
ence of Bd mesons, which is predicted to be negligible in the SM and many of its CPT 
conserving extensions, may reveal indirect CPT violation. Time dependence of the di- 
rect CP violating asymmetry for flavour specific decays, which is time independent and 
vanishes for decays dominated by only one amplitude may establish CPT violation as 
well as a large life time difference. The theoretically clean inclusive semileptonic decays 
having relatively large branching ratios might be particularly suitable in this context. 
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Figure 1: The time evolution of the time dependent asymmetry (acp(t)) f° r I m <5 = 0.1 
for any one amplitude dominated process (in all figures t is in units of rg/lO, where Tb is 
the average B° life time ). See text for the details. 
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Figure 2: The variation of Dft) = ProhfR — > f) - ProhfR — > f) as function of time for 
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Figure 3: The variation of S(t)= Prob(B — > f) + Prob(B — > f) as function of time. 
This quantity is sensitive to Re 5 only. For comparison we have plotted for Ke5=0.0 (SM) 
(the dotted curve) and Re<5 = 0.5 (the solid curve). 
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Fiemre 4: The variation of the time integrated asvmmetrv a for a one amnlitude 



